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We calculate the momentum distribution and the second-order correlation function in momentum 
space, g^ (p, p', t) for molecular dimers that are coherently formed from an ultracold atomic gas 
by photoassociation or a Feshbach resonance. We investigate using perturbation theory how the 
quantum statistics of the molecules depend on the initial state of the atoms by considering three 
different initial states: a Bose-Einstein condensate (BEC), a normal Fermi gas of ultra-cold atoms, 
and a BCS-type superfluid Fermi gas. The cases of strong and weak coupling to the molecular field 
are discussed. It is found that BEC and BCS states give rise to an essentially coherent molecular 
field with a momentum distribution determined by the zero-point motion in the confining potential. 
On the other hand, a normal Fermi gas and the unpaired atoms in the BCS state give rise to a 
molecular field with a broad momentum distribution and thermal number statistics. It is shown 
that the first-order correlations of the molecules can be used to measure second-order correlations 
of the initial atomic state. 
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I. INTRODUCTION 

The basis of some of the most exciting developments 
in ultra-cold atomic physics in recent years has been 
the use of Feshbach resonances 0, and photoasso- 
ciation 0, H to tune the strength of the interactions 
between atoms as wells as to create ultracold diatomic 
molecules startingfrom an ultracold atomic gas of bosons 
or fermions 0, MM 0, H H E3 The availability of tun- 
able interactions has made possible the study of many 
model systems of condensed matter theory in a very 
controlled fashion [H El El El El E3 In particu- 
lar, the BEC-BCS crossover, which predicts a continuous 
transition from a BCS superfluid of atomic fermions to 
a molecular BEC as the interaction strength is varied 
from attractive to repulsive, has attracted a considerable 
amount of attention, both experimentally and theoreti- 
cally [S El El EldS E3 

A difficulty in the studies of the BEC-BCS crossover 
has been that they necessitate the measurement of higher 
order correlations of the atomic system. While the mo- 
mentum distribution of a gas of bosons provides a clear 
signature of the presence of a Bose-Einstein condensate, 
the Cooper pairing between fermionic atoms in a BCS 
state hardly changes the momentum distribution or spa- 
tial profile as compared to a normal Fermi gas. This 
poses a significant experimental challenge, since the pri- 
mary techniques for probing the state of an ultracold gas 
are either optical absorption or phase contrast imaging, 
which directly measure the spatial density or momentum 
distribution following ballistic expansion of the gas. In 
the strongly interacting regime very close to the Fesh- 
bach resonance, evidence for fermionic superfluidity was 
obtained by projecting the atom pairs onto a molecular 



state by a rapid sweep through the resonance [rH ll^| . 
More direct evidence of the gap in the excitation spectra 
due to pairing was obtained by rf spectroscopy [22] and 
by measurements of the collective excitation frequencies 

Still, the detection of fermionic superfluidity in the 
weakly interacting BCS regime remains a challenge. The 
direct detection of Cooper pairing requires the measure- 
ment of second-order or higher atomic correlation func- 
tions. Several researchers have proposed and imple- 
mented schemes that allow one to measure higher order 
correlations |l7l l25l 12a, 1271 12a . |29| but those methods are 
still very difficult to realize experimentally. 

While the measurement of higher order correlations is 
challenging already for bosons, the theory of these corre- 
lations has been established a long time ago by Glauber 
for photons I30L l3ll l32 | . For fermions however, despite 
some efforts [33 a satisfactory coherence theory is still 
missing. 

To circumvent these difficulties we suggested in an ear- 
lier publication [3J|, guided by the analogy with three- 
wave mixing in classical optics, to make use of the non- 
linear coupling of atoms to a molecular field by means of 
a two-photon Raman transition or a Feshbach resonance. 
The nonlinearity of the coupling links first-order correla- 
tions of the molecules to second-order correlations of the 
atoms. Furthermore the molecules are always bosonic so 
that the well-known coherence theory for bosonic fields 
can be used to characterize them. Considering a simpli- 
fied model with only one molecular mode, it was found 
that the molecules created that way can indeed be used 
as a diagnostic tool for second-order correlations of the 
original atomic field. Naturally, due to the restriction to 
a single mode, the information one can gain about the 
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atomic state is very limited. 

In this paper we extend the previous model to take 
into account all modes of the molecular field, the hope 
being that in doing so, more detailed information about 
the atomic state can be obtained. Specifically, we cal- 
culate the momentum distribution of the molecules and 
the normalized second-order correlation function, g^ for 
different momentum states of the molecules using pertur- 
bation theory. We consider the limiting cases of strong or 
weak atom-molecule coupling as compared to the relevant 
atomic energies. The molecule formation from a Bosc- 
Einstein condensate (BEC) serves as a reference system. 
There we can rather easily study the contributions to the 
molecular signal from the condensed fraction as well as 
from thermal and quantum fluctuations above the con- 
densate. The cases of a normal Fermi gas and a BCS 
superfluid Fermi system are then compared with it. We 
show that the molecule formation from a normal Fermi 
gas and from the unpaired fraction of atoms in a BCS 
state has very similar properties to those of the molecules 
formed from the non-condensed atoms in the BEC case. 
The state of the molecular field formed from the pairing 
field in the BCS state on the other hand is similar to that 
resulting from the condensed fraction in the BEC case. 
The qualitative information gained by the analogies with 
the BEC case help us gain a physical understanding of 
the molecule formation in the BCS case where direct cal- 
culations are difficult and not nearly as transparent. 

This paper is organized as follows: In section [n] we 
introduce the model Hamiltonian used to describe the 
coupled atom-molecule system. In sections IIIII to [V] we 
present the calculations of momentum distribution and 
second factorial moment of the molecular field for a BEC, 
a normal Fermi gas and a BCS-type state, respectively, 
considering the cases of strong and weak coupling in each 



Here £k = k 2 /2M is the kinetic energy of an atom of 
mass M and momentum k and = + v is the 
energy of a molecule with momentum k and detuning 
parameter v. c^ a and cj^ are fermionic annihilation and 
creation operators for plain waves in quantization vol- 
ume V with spin a. t/ tr (k) = V~ 1 ^ 2 J v d 3 xe~ lkx U t r(x) 
is the Fourier transform of the trapping potential Uti{v) 
and Uq = Aita/M is the background scattering strength, 
g is the effective coupling constant of the atoms to the 
molecules and we use units with % = 1 throughout. We 
assume that the trapping potential and background scat- 



case. Details of the calculations are given in the appen- 
dices El and H2 



II. MODEL 

The general procedure that we have in mind is the fol- 
lowing: The atomic sample is prepared in some initial 
state, whose higher order correlations we seek to ana- 
lyze. At time t = the coupling to a molecular field is 
switched on. While the initial atomic state corresponds 
to a trapped gas, we assume that the molecules can be 
treated as free particles. This is justified if the atomic 
trapping potential does not affect the molecules, or if 
the interaction time between the atoms and molecules is 
much less than the oscillation period in the trap. Finally, 
the state of the molecular field is analyzed by standard 
techniques, e.g. time of flight measurements. 

We consider the three cases where the atoms are 
bosonic and initially form a BEC, or consist of two species 
of ultra-cold fermions (labeled by u =T, |), with or with- 
out superfluid component. In the following we describe 
explicitly the situation for fermions, the bosonic case be- 
ing obtained from it by omitting the spin indices and 
by replacing the Fermi field operators by bosonic field 
operators. 

Since we are primarily interested in how much can be 
learned about the second-order correlations of the initial 
atomic cloud from the final molecular state, we keep the 
physics of the atoms themselves as well as the coupling 
to the molecular field as simple as possible. The coupled 
fcrmion-molecule system can be described by the Hamil- 
tonian dium 



(1) 



tering are relevant only for the preparation of the initial 
state before the coupling to the molecules is switched on 
at t = and can be neglected in the calculation of the 
dynamics. This is justified if g^fN 3> Uon,uJi where n is 
the atomic density, N the number of atoms, and are 
the frequencies of Ut r (r) that is assumed to be harmonic. 
In experiments, the interaction between the atoms can 
effectively be switched off by ramping the magnetic field 
to a position where the scattering length is zero. 

Regarding the strength of the coupling constant g, two 
cases are possible: g^/N can be much larger or much 



k,<r k ki ,k2 ,cr 

+ ^ H C+qT^L-qi^icidT +9 XX 5 i/ 2 + k i fi q/2-kT + H - C - 

q,ki,k 2 \ q,k 



r 



3 



smaller than the characteristic kinetic energies involved. 
For fermions the terms broad and narrow resonance have 
been coined for the two cases, respectively, and we will 
use these for bosons as well. Both situations can be real- 
ized experimentally, and they give rise to different effects. 
We examine both limiting cases and use the suggestive 
notation E^ n <C g\/N and Ekin ^> gVN for the two 
cases, where -Ekin denotes the characteristic kinetic en- 
ergy of the atoms. It corresponds to zero point motion 
for condensate atoms, to the thermal energy, k B T for 
non-condensed thermal bosons, and to the Fermi energy 
for a degenerate Fermi gas. 

Our analysis is based on the assumption that first order 
time-dependent perturbation theory is applicable. This 
requires that the state of the atoms does not change sig- 
nificantly and consequently, only a small fraction of the 
atoms are converted into molecules. It is reasonable to 
assume that this is true for short interaction times or 
weak enough coupling. Apart from making the system 
tractable by analytic methods there is also a deeper rea- 
son why the coupling should be weak: Since we ulti- 
mately wish to get information about the atomic state, 
it should not be modified too much by the measurement 
itself, i.e. the coupling to the molecular field. Our treat- 
ment therefore follows the same spirit as Glauber's origi- 
nal theory of photon detection, where it is assumed that 
the light-matter coupling is weak enough that the detec- 
tor photocurrent can be calculated using Fermi's Golden 
rule. 



is the condensate wave function in the Thomas-Fermi 
approximation and c is the annihilation operator for an 
atom in the condensate, Rtf — (15A r a/a osc ) 1 / 5 a osc is the 
Thomas-Fermi radius, N the number of atoms, a is their 
scattering length and a osc is the oscillator length of the 
atoms in the trap. In accordance with the assumption of 
low temperatures and weak interactions we do not distin- 
guish between the total number of atoms and the number 
of atoms in the condensate. The fluctuations S"ijj{x) are 
small and those with wavelengths much less than Rtf 
will be treated in the local density approximation while 
those with wavelengths comparable to Rtf can be ne- 
glected, mm nt 



A. Broad resonance, i?km <C g^/N 



We are interested in the momentum distribution of the 
molecules 



n(p,t) = (aUt)a p (t)) 



(4) 



III. 



BEC 



which for short times, t, can be calculated using per- 
turbation theory: We expand n(p, t) in a Taylor series 
around t = and make use use of the Heisenberg equa- 
tions of motion 



We consider first the case where the initial atomic state 
is a BEC in a spherically symmetric harmonic trap. We 
note that all of our results can readily be extended to 
anisotropic traps by an appropriate rescaling of the coor- 
dinates in the direction of the trap axes. We assume that 
the temperature is well below the BEC transition tem- 
perature and that the interactions between the atoms are 
not too strong. Then the atomic system is described by 
the field operator 



^(x) =Xo(x)c + ^(x), 



where 



Xo(x) 
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. da p (t) 

1 dt 
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and similarly for a p . Here we have neglected the kinetic 
energy term, a step that is legitimate for a broad reso- 
nance since the interaction energy gVN is much larger 
than the difference in the kinetic energies between the 
atoms and molecules. Consequently for short enough 
times, t < (gy/N)^ 1 , energy conservation can be vio- 
lated in the formation of molecules in a fashion similar 
to the Raman-Nath regime of atomic diffraction. 

To lowest non- vanishing order in gt wc find 
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™BEC,b(p,t) = (gt) 2 \ 6 p/2-k 1 6 p/2+k 1 £ p/2+k 2 Cp/2-k 2 ) +0((gt) 4 ), (6) 



where the atomic operators are the initial (t = 0) oper- 



ators. This expression can be evaluated by making use 
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of the decomposition of the atomic field operator and 
the local density approximation for the part describing 
the fluctuations. The details of this calculation are given 



J 



in appendix ^ To first non- vanishing order in the fluc- 
tuations wc find 



n BEC , h (p,t) = (gt) 2 N(N-l)V\x 2 (p)\ 2 + (gt) 2 4N J ^ (x)«5c p (x)) 



(7) 



From this expression we see that our approach is justified 
if (VN gt) 2 <C f because for such times the initial atomic 
state can be assumed to remain undepleted. 

In the local density approximation the expectation 
value (<5c p (x)(Sc p (x)) for the number of fluctuations with 
momentum p at x can be evaluated by assuming that at 
each x we have a homogenous BEC with density n(x) and 
using the Bogoliubov transformation to quasi-particle op- 
erators, dk(x). One finds 

((Set (x)(JCp(x)) = W 2( x ) + ( w 2 (x)+w 2 (x)) (x)dp(x)) , 

(8) 

(9) 
(10) 

(11) 



with Bogoliubov amplitudes 
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^p(x) 
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e P (x) 



and quasi-particle energies 



e p (x) = Jet + 2e p n(x)C/ . 



The quasi-particle distribution is given by a thermal Bose 
distribution 



(dp(x)a p (x)) = 



1 



e e~ p (x)/fc B T _ I ■ 



(12) 



The momentum distribution (jTJ is illustrated in Fig. 
H The contribution from the condensate is a collective 
effect, as indicated by its quadratic scaling with the atom 
number. It clearly dominates over the incoherent contri- 
bution from the fluctuations, which is proportional to the 
number of atoms. The momentum width of the contri- 
bution from the condensate is roughly 2tt/ R^f which is 
much narrower than the contribution from the fluctua- 
tions, whose momentum distribution has a typical width 
of l/£, where £ = (Snan)" 1 / 2 is the healing length. Us- 
ing the Thomas-Fermi wave function for the condensate, 
Eq. ©, we can calculate the condensate contribution in 
closed form as 



225N(N~l)(gt) 2 

nBEC ' b(p ' i) = ^R^f 



6 sin pRtf 6 cos pRtf 
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The terms of order N are corrections due to the non- 
condensed part. 

Using the same approximation scheme we can calculate 
the second-order correlation 



g {2) (pi,t 1 ;p 2 ,t 2 ) 



(a P1 fti)a p2 ft 2 )a P2 (<2)a Pl (*i) 
n(p 1 ,t 1 )n(p 2 ,t 2 ) 



If we neglect fluctuations we find 



(14) 
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(iV-2)! 2 
(N-A)\Nl 

l-±+0(N- 2 ). (15) 



For iV — > oo this is very close to 1, which is characteristic 
of a coherent state. This result implies that the number 
fluctuations of the molecules are very nearly Poissonian. 
The fluctuations lead to a larger value of g^ 2 \ making 
the molecular field partially coherent, but their effect is 
only of order 0{N^ 1 ). 

The physical reason why the resulting molecular field 
is almost coherent is of course clear: The condensed frac- 
tion of the atomic field operator is dominant. In expecta- 
tion values the operators c and & take on values \/N — n, 
with a number n <C N depending on the position of the 
operator in the expectation value. When dividing by the 
normalizing expectation values n(p, f), \/N — n can be 
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FIG. 1: (color online) Momentum distribution of molecules 
formed from a BEC (red dashed line) with a — Clow 
and T — 0.1T C and a BCS type state with fed = 0.5 and 
(Jose = 5fc^ 1 (0) (blue solid line), both for TV = 10 5 atoms. 
The BCS curve has been scaled up by a factor of 20 for easier 
comparison. The inset shows the noise contribution for BEC 
(red dashed) and BCS (blue) case. The latter is simply the 
momentum distribution of molecules formed from a normal 
Fermi gas. The local density approximation treatment of the 
noise contribution in the BEC case is not valid for momenta 
smaller than 2n/£ (indicated by the red dotted line in the 
inset). Note that the coherent contribution is larger than the 
noise contribution by five orders of magnitude in the BEC 
case and three orders of magnitude in the BCS case. 

replaced by y/W with accuracy 0{N~ X ) and hence c and 
can be replaced by yf~N independent of their position 
in the expectation value. The field operator can thus be 
replaced by a c- number field ^>(x) — > y/Nxo(x), the mean 
field, which explains the almost perfect factorization of 
the correlation functions. 

Another way to understand this is to consider the 
single-mode BEC state, |*(0)) = (c^W 2 ) |Q)/VM. 
The coupling to the molecular field will transform this 
state into w (aa) + 0&)( N / 2 ) |0}/VM, where 

a <C 1. This leads to a binomial distribution for the 
number of molecules. In the limit that N — > oo the Bi- 
nomial distribution goes over to the Poisson distribution. 



B. Narrow resonance, Skin 3> gy/~N 

For a narrow resonance, the typical kinetic energies as- 
sociated with the atoms and molecules, -Ekin, are much 
larger than the atom-molecule interaction energy. This 
implies that even for very short interaction times, t < 
(gy/N)" 1 , the phase of the atoms and molecules can 
evolve significantly, E\^ n t 3> 1. Consequently, only tran- 
sitions between atom pairs and molecules that conserve 
energy can occur. In this case, it is convenient to go over 
to the interaction representation 

Cp(t) -» e- ie "'c p (i), o p (i) e- iE ^a p (t), (16) 
where for notational convenience, we will denote the in- 
teraction picture operators by the same symbols as the 
Heisenberg operators used in the previous subsection. 

The equations of motion in the interaction picture, 

i^M = g ^ e^K-p/— p/=-)*c p/2+k c p/2 _ k , 

(17) 

can be approximately integrated by treating the atomic 
operators as constants, leading to 

%>(*) = 2J ^( E P ~ £ p/2+k - Cp/2-k, £)Cp/2+kCp/2-k7 

(18) 

where we have introduced 

iujt _ ^ 

A(w,t)=0lim— — . (19) 

rjj.0 iuj + rj 




The condition under which this step is justified is ana- 
lyzed below. As in the broad resonance case we can insert 
this expression in n(p,i). The calculation of the result- 
ing integrals over expectation values of the atomic state 
is however considerably subtler than in for the broad res- 
onance case and is presented in details in appendix^] In 
the limit !/(> 1 wc find 



»BBC,n(p,*) = N(N-1) 



V 2 M 3 g 2 
16tt 4 
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Jo V v- 



x J dzxo (^\Jp 2 /4: + Mui - py/TEjz^ xo (\/ P 2 / 4 + Mui + ps/Wu 



3g 2 tV^APR 3 TF [fx ^ 

iyo P /2,vrn; 87r2 



y (5c p (x)«5c p (x)). 



(20) 



r 



In the second term in (|20|l we have defined 



As before, the contribution from the condensate is clearly 



4ttV 



3Rj^p j —i 



dz\xo(Vp 2 + P' 2 ~ 2pp'z) j 



0(1), \p~ P '\<2tt/R tf 



(21) 
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FIG. 2: In the narrow resonance case for momenta p much 
larger than the momentum width of the condensate one atom 
is taken out of the condensate and the other atom is taken 
from the non-condensed part and has momentum close to p 
in order for momentum conservation to be satisfied. Since 
the total energy of the atoms has to match the total energy 
of the molecule, only molecules with momenta 2 V Mv can be 
formed for each detuning v. 




FIG. 3: Momentum distribution of molecules formed from a 
BEC of N = 10 5 atoms with scattering length a = 0.01a OS c 
at T = 0.1T C for a narrow resonance. 



dominant. The integral in the first term in Eq. (|20() is 
proportional to the amplitude for finding an atom pair 
with center of mass momentum p and total kinetic energy 
v. Because xo drops to zero on a scale of 2it/Rtf this 
amplitude is essentially zero if p > 27t/_Rtf or v > M n Rii . 

The second term in Eq. I|20() originates from molecules 
that are formed from an atom in the condensate and 
a non-condensed atom. Since the atom momentum 
< 27t/_Rtf in the condensate is very small compared 
to the momentum |p| ~ l/£ of a non-condensed atom, 
the molecular momentum is essentially due to the non- 
condensed atom. On the other hand, energy conservation 
implies that v+p 2 /AM » p 2 /2M if p > 2ir/R TF . Conse- 
quently for a given detuning v, molecules with momenta 
in a shell of radius 2\J Mv and width 2tt/Rtf are formed 
from one atom in the condensate and another atom taken 
from the non-condensed part with a momentum that also 
lies in a spherical shell in momentum space around p with 
thickness 2tt/Rtf- Momentum and energy conservation 
are illustrated in Fig. [3 Figure |3| shows a typical exam- 
ple for the momentum distribution. 

Equation (|20|l allows us to extract the criterion for the 
applicability of our approximation scheme, i.e of treat- 
ing the atomic state as being undepleted. The coherent 
contribution will only be nonzero if |p| < 2tt/Rtf and 
for these momenta they can be neglected, as we have 
seen. Requiring that the number of molecules remains 
much smaller then the initial number of atoms leads to 
the condition 



1 



(22) 



In the opposite case |p| > 2tt/Rtf the coherent contri- 
bution is essentially zero and we need only consider the 
incoherent contribution. Requiring that the number of 
molecules with momentum p be much smaller than the 



number of non-condensed atoms with that same momen- 
tum leads to 



gt < N~ 



1 



(23) 



IV. NORMAL FERMI GAS 

For a normal Fermi gas we restrict ourselves to the 
case of zero temperature, T = 0. For temperatures T 
well below the Fermi temperature Tp, the corrections to 
our results are of order (T/Tf) 2 or higher, and do not 
lead to any qualitatively new effects. 

Again, we treat the gas in the local density approxi- 
mation where the atoms locally fill a momentum sea 



\NFG) = J] 



4lo) 



(24) 



|k|<fc F (x) 



with local Fermi momentum fcp(x) and |0) being the 
atomic vacuum. The Fermi momentum is related to the 
local chemical potential 



A*ioc(x) = fx - i/tr(x) 



by means of 



Mloc(x) 



„(x) _ (3^ 2 n(x)) 2 / 3 



2M 



2M 



(25) 



(26) 



Here, (j,q = (37r 2 no) 2 / 3 /(2M) is the chemical potential of 
the trapped gas, and no is the density at the center of 
the trap for each of the spin states. The Hartree-Fock 
mean field has been neglected because it gives rise only 
to minor corrections and doesn't lead to a qualitatively 
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new behavior. The density distribution of the trapped 
gas is given by the Thomas-Fermi result [40| . 



where we have introduced the local density of atom pairs 
with center of mass momentum p, 



N 8 



-, 3/2 



R 2 F 



(27) 



(487V) 



1/6, 



is the Thomas-Fermi radius 



where Rf = 
for fermions. 

Using the same perturbation methods as described in 
the previous section for bosons, we can calculate the mo- 
mentum distribution of the molecules and their correla- 
tion function </ 2 ' for a broad and for a narrow resonance 
by first calculating the density of the desired quantity 
at a position x and then integrating the result over the 
volume of the gas. 



A. Broad resonance E^ ln <C gVN 

To deal with the case of fermions, we modify Eq. © 
by reintroducing the spin of the atoms. The integral over 
the relative momentum of the atom pairs can be carried 
out exactly to give 



"NFG,b(P,X, t) 



( ff t) 2 ^(p,x), |p|<2fc F (x) 



0. 



|p| > 2k F (x), 



(28) 



*UP,x) 



nk% (x) 
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k F (x) 



k F (x) 



(29) 

-Ffc(p,x) can be visualized as the integration over the in- 
tersection of two Fermi seas shifted by p relative to each 
other, as depicted in Fig. 01a). 

The characteristic width of the momentum distribu- 

1/3 

tion of the molecules is kp oc n Q which is typically 
much wider than the distribution found in the BEC case. 
From Eq. (|28|l the number of molecules produced scales 
linearly with the number of atoms. This is because in 
contrast to the BEC case, the molecule production is a 
non-collective effect. Each atom pair is converted into a 
molecule independently of all the others and there is no 
collective enhancement. The integration of these results 
over the volume of the cloud can easily be done numeri- 
cally and is shown in the inset in Fig. ^ 

Similarly, we can calculate the local value of at 
position x, from Eq. (|14fl . 



5ioc(Pi^i;P2^2) = 



S.F b ( Pl) x)F b (p 2 ,x)- J <2krc(p 2 /2 + k)rc(p 2 /2-k)n( Pl - p 2 /2 - k) 
J dkn(p 2 + k - pi/2)ra(pi/2 + k)n(pi/2 - k) 

J dkn(p 2 - Pl /2 + k)n(pi/2 - ki)n(p 2 /2 - k 2)j/ F 6 ( Pl , x)F 6 (p 2 , x). 



(30) 



This result simplifies considerably for Pi = P2 = P; 
0ioc(P> x >*) ee ggl(p,t;p,t,x) 

= 2 ( 1 -^—] . (31) 

As in the case of the BEC, the time dependence in 
g( 2 ) (p, x, t) cancels at this level of approximation. How- 
ever, in contrast to the case of a BEC there is some de- 
pendence on the momentum left. 

The origin of the factor of two in g^> (p, x, t) is the 
following: The two molecules that are being detected in 
the measurement of g^ can be formed from four atoms 
in two different ways and the two possibilities both give 
the same contribution. Eq. I|31|l indicates that the statis- 
tics of the molecules are super-Poissonian, similarly to a 
thermal field. 

By the following argument we can convince ourselves 
that not only the second-order correlations look thermal, 



but that the entire counting statistics of each momentum 
mode is thermal. Each molecular mode characterized 
by the momentum p is coupled to a particular subset 
of atom pairs selected by momentum conservation. In 
the short time limit, each atom pair with center of mass 
momentum p is converted into a molecule in the corre- 
sponding molecular mode independently of all the other 
atom pairs and with uncorrelated phases. Thus we ex- 
pect the number statistics of each molecular mode to be 
similar to that of a light field in thermal equilibrium with 
a reservoir with which there is an incoherent exchange of 
energy. 

B. Narrow resonance, Euin 3> ffV^V 

In this case, the molecules formed have to satisfy en- 
ergy and momentum conservation, as illustrated in Fig. 
0Jb). We are then lead to a calculation very similar to 



8 




the one presented in appendix[B]for the contribution from 
the non-condensed fraction of atoms for the BEC case. 
After integrating over the volume of the cloud we find 




FIG. 4: (a) Illustration of the number of atom pairs with 
center of mass momentum q. These atoms can be transformed 
into a molecule of momentum q if energy conservation plays no 
role, (b) Density of atom pairs with center of mass momentum 
q and total kinetic energy 2V Mv. These atom pairs can be 
converted into a molecule with momentum q and detuning v. 



The number of molecules produced is proportional to the are the same as those in the broad resonance case except 

number of atom pairs that satisfy momentum and energy that the integration is now over pairs of atoms that also 

conservation and hence scales linearly with the number satisfy energy conservation. For the particular case of 

of atoms, indicating that molecule formation is not a col- Pi = P2 = P> we obtain an expression with the exact 

lective effect. Figure0shows the momentum distribution same form as Eq. ()31|) except that F(,(p,x) must be 

for typical parameters. It is much wider than the momen- replaced with F n (p, v, x), 
turn distribution for the BEC case in both momentum 
space and in energy width. 

We don't give the lengthy and complicated expression 
for fl£,c(pi,ii,P2)t2) because its qualitative properties 




(32) 



F n (p,v, x) 




(33) 




V. 



BCS STATE 



Let us now consider a system of Fermions with attrac- 
tive interactions, Uq < 0, at temperatures well below the 
BCS critical temperature. As is well known, for these 
temperatures the attractive interactions give rise to cor- 
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FIG. 5: Momentum distribution of molecules produced from a 
normal Fermi gas in the narrow resonance limit for g = 10 _3 /x. 



Here £k(0 = £k — Mioc(^) is the kinetic energy of an atom 
measured from the local chemical potential defined as 



Hioc{r) = Mo - U(r) - U n(r). 



(37) 



In contrast to the normal Fermi gas, we have included 
a Hartree-Fock mean-field energy to the local chemical 
potential since we can no longer ignore the effect of the 
two-body interactions in the gas. To an excellent approx- 
imation we can use the relation Eq. (|26|l between density 
and local chemical potential. Then, for a given number 
of atoms TV, Eq. (gZj) is an implicit equation for /iq. We 
solve it numerically and hence determine the density pro- 
file n(r) and the local chemical potential fj,i oc (r). 



relations between pairs of atoms in time reversed states 
known as Cooper pairs. We assume that the spherically 
symmetric trapping potential is sufficiently slowly vary- 
ing that the gas can be treated in the local density ap- 
proximation. More quantitatively, the local density ap- 
proximation is valid if the size of the Cooper pairs, given 
by the correlation length 

A(r) = VF(r)/TtA(r), 

is much smaller than the oscillator length for the trap. 
Here, Vf{t) is the velocity of atoms at the Fermi surface 
and A(r) is the pairing field at distance r from the origin, 
which we take at the center of the trap. 

Before turning to the coupled atom-molecule system 
we outline our treatment of the atomic system. We 
closely follow the approach of Houbiers et. al. |4l|. We 
assume that locally at each r, the wave function can be 
approximated by the BCS wave function for a homoge- 
nous gas, 



\BCS(r)) = ]J(u k (r)+v k (r)cl_ 



with Bogoliubov amplitudes 



k,T C ka 



)|0>, (34) 



The gap parameter A(r) = Uq/2V J2k u k{ r ) v it( r ) i s 
determined by the gap equation 



2fc F (0)a 



= Mo^ 3 (0) 



dkk 2 



1 



V£(r)+A2(r) 



&(r) 
(38) ' 



where the ultra-violet divergence has been removed by 
renormalizing the bare background scattering strength 
to the two-body T-matrix using the Lippmann-Schwinger 
equation (see ref. |41|). We solve the gap equation nu- 
merically using the previously determined local chemical 

A*loc- 



A. Broad resonance, Ekin <C gVN 




&(r) 



(35) We find the momentum distribution of the molecules 
from the BCS type state by repeating the calculation 
done in the case of a normal Fermi gas. For the BCS wave 

(36) function, the relevant atomic expectation values factorizc 
as 



C p/2-k 1 ,T C p/2+k 1 ,l C P/ 2 + k 2^T C p/2-k 2 a 



S P /2- kl ,T S p/2+k x ,J. / \ c p/2+k 2 ,lC p /2-k 2 ,T 
C p/2-k 1 ,T C P/ 2 +k2.l) \ C P /2+k 1 a C P/2-k 2 ,T 



(39) 
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and the momentum distribution of the molecules be- 
comes 



J 



n.Bcs,b(p,i) = (gt) 2 
« (gt) 2 



E 



-p/2+k,l C p/2-k : T 



k 



+ E( £ P/2- l 
k 

+ "NFG,b(P,i)- 
I 



M Cp/2+k,l ) \C p/2+ka Cp/2-k,T 



(40) 



In going from the first to the second line we have as- 
sumed that the interactions are weak enough so that the 
momentum distribution of the atoms is essentially that 
of a two component Fermi gas. This is justified because 
the Cooper pairing only affects the momentum distribu- 
tion in a small shell of thickness 1/X(r) <C kp(r) around 
the Fermi surface. 



J 



The first term involves the square of the pairing field. 
It is proportional to the square of the number of paired 
atoms which, below the critical temperature, is a finite 
fraction of the total number of atoms. This quadratic 
dependence indicates that it is the result of a collective 
effect. This term can be related to the two-point corre- 
lation function in position space as 



c p/2+k,lCp/2-k,T 



d 3 xd 3 i 
V 



•<Vi(x-r/2)Vi T (x + r/2)) 



(41) 



r 



where ifo ± are the atomic field operators in position 
space. In the local density approximation, the correla- 
tion function varies with x on a length scale -Rtf- On 
the other hand, the expectation value in the integral falls 
off to zero for r > A(r) and we can therefore treat the 
correlation function as being independent of x when per- 
forming the integration over r, 



'p/2+k,J.Cp/2-k,T 



d 3 x 



-ipx 



Ck,l c -kT/ 



(42) 



The expectation value on the right hand side is evalu- 
ated using the local density approximation at position x. 
Inserting the result into Eq. (|40|1 and making use of the 
gap equation we find 



«BCS,b(p,t) = (gtf 



d 3 xe~ l * p 1 



2aA\ 2A(z) 



r 



^NFG,b(P, t) 



(43) 



Following ref. [4l| we have replaced the bare background 
coupling strength by 



U = U - 



1 



2aA 



(44) 



where A is a momentum cut-off and is of the order of the 
inverse of the range of the inter-atomic potential. 



Using the numerically determined A(x) we can read- 
ily perform the remaining Fourier transform in Eq. I|43|) . 
The result of such a calculation is shown in Fig. Since 
the gap parameter changes over distances of order Rtf 
the contribution from the pairing field has a typical width 
of order 1/Rtf- This is very similar to the BEC case. 
The background from the unpaired atoms on the other 
hand has a typical width kp (0) = (37r 2 7io) 1 / 3 which is 
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similar to the width of the noise contribution in the BEC 
case, see inset in Fig. Q] This similarity can be under- 
stood by recalling that for a weakly interacting conden- 

1 /3 

sate, n a = (3 <C 1, which when substituted into the 
definition of the healing length gives l/£ = (87r/3) 1 / 2 nQ . 
For typical j3 ~ 0.1, this is comparable to fcp(0) for equal 
densities. 

Because of the collective nature of the coherent con- 
tribution it will dominate over the background, riNFG,b, 
for strong enough interactions and large enough particle 
numbers. The narrow width and the collective enhance- 
ment of the molecule production are the reasons why the 
momentum distribution of the molecules is such an excel- 
lent indicator of the presence of a superfluid component 
and the off-diagonal long range order accompanying it. 

For weak interactions such that the coherent contribu- 
tion is small compared to the incoherent contribution, the 
second order correlations are close to those of a normal 
Fermi gas given by Eq. I|31|) . g^(p, x, t) « 2. However, 
in the strongly interacting regime, kp\a\ ~ 1, and large 
N , the coherent contribution from the paired atoms dom- 



inates over the incoherent contribution from unpaired 
atoms. In this limit one finds that the second-order cor- 
relation is close to that of the BEC, (p, x, t) « 1. The 
physical reason for this is that at the level of evenorder 
correlations the pairing field behaves just like the mean 
field of the condensate. This is clear from the factor- 
ization property of the atomic correlation functions, Eq. 
(|39|l . in terms of the normal component of the density 
and the anomalous density contribution due to the mean 
field. In this case, the leading order terms in N are given 
by the anomalous averages. In the strongly interacting 
limit, the contribution from the 'unpaired' atoms is very 
similar in nature to the contribution from the fluctua- 
tions in the BEC case. 



B. Narrow resonance, i5km 3> g\/~N 

A calculation similar to the one presented in Appendix 
|5]for the BEC case leads to 



nBCS,n(p,£) = 



J2 A (V~ fc 2 / M )(Cp/2+k,l 5 p/2-k,T/ 



' ^NFG,n(P,i), 



(45) 



where we have assumed again that the gas is weakly in- and performing similar manipulations as in the broad 
teracting. Inserting Eq. I|19|l for A in the limit vt — > oo resonance case leads to 



nBCS,n(P,*) 



g 2 M 3 



1 2 



,oo/ 1 

/ duj'Juj n5(v — u>) + iP 

Jo V v- 



drr sm(pr) (c^ | c_ y/m< f ) 



+ "-NFG.n 



(46) 



where again the pairing field (ck t iC-k,i)\ r = Uk(r)vk(r) 
can be evaluated using the local density approximation. 
Figure H3 shows an example of the pairing field across the 
trap. Two qualitatively different cases have to be distin- 
guished depending on the strength of the interactions. 

If the interactions are fairly strong so that the pair- 
ing field, Uk{r)vk{r), is nonzero in a rather wide region 
around kp(r), the pairing field will be a slowly varying 
function across the atomic cloud. Then the remaining 
integral in eq. 1)46(1 can be easily evaluated numerically 
and we find a momentum distribution of the molecules 
which is similar to the BEC case. This limit is illustrated 



in fig. \7\ The width of the momentum distribution of 
the molecules is again of order 1/Rtf- It is known that 
in the strongly interacting limit, the size of the Cooper 
pairs becomes comparable to the interparticle spacing, 
A(r) ~ l/fcj?(r). The Cooper pairs are no longer delo- 
calized across the extent of the cloud but now approach 
the limit of localized bosonic " quasi-molecules" . Thus 
it is not surprising that the momentum distribution of 
molecules formed from a BCS type state approaches the 
one we found in the BEC case. 

On the other hand, if the interactions between the 
atoms are weak the pairing field is a very narrow func- 
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k [k (0)j 



30 r[1/L(0)] 50 



FIG. 6: Pairing field (ck,|C_k,|) = utVk across the trap for 
kFd = 0.5 and a osc = 5k F 1 (0). The black solid line indicates 
the local Fermi momentum kp(r). 



tion of momentum and hence, for fixed momentum, also 
of position as can be seen from Fig. [fJl Then the integral 
in Eq. I|46l) has contributions from a rather narrow re- 
gion in space only and the momentum distribution will 
accordingly become wider and smaller. 

Probing the BCS system in the narrow resonance 
regime also yields spatial information about the atomic 
state. By tuning v = 2fi\ oc (r) the molecular signal is 
most sensitive to the pairing field near r and less sensi- 
tive to other regions in the trap. 

A qualitative difference between the BCS and BEC 
cases becomes apparent if one looks at the number of 
molecules as a function of the detuning. While we find 
that there is only a very narrow distribution of detunings 

with width ~ < - 27r ^J F - ) that leads to molecule formation 
in the BEC case, molecules are being formed for detun- 
ings well below v ~ 2/io- The non- homogeneity of the 
trapped atom system manifests itself in a completely dif- 
ferent way in the two cases. In the BEC case, the total 
energy of a particle in the condensate is just the chemi- 
cal potential while the kinetic energy of a particle is very 
small compared to the mean field energy in the Thomas 
Fermi limit. Upon release, the atoms in the condensate 
all have a spread in kinetic energies that is of the or- 
der ^ 2lr ^ f ^ due entirely to zero point motion. On the 
other hand, in the BCS case the superfluid forms at each 
position near the local Fermi momentum. Hence, atoms 
in the BCS state have a large energy spread that is of 
order ~ /iq with the kinetic energies of the paired atoms 
being centered around ^io C (?") with a width A(r). 

For the second order moment, g^ (p, x, t), the same 
general arguments that were put forward in the discus- 
sion of the broad resonance also apply to the narrow res- 
onance. In the strongly interacting limit, the molecular 
field is again approximately coherent with a noise contri- 
bution from the unpaired fermions. 



FIG. 7: Momentum distribution of molecules formed from a 
strong coupling (fc^a = 0.5) BCS system in the narrow reso- 
nance regime as a function of detuning v for a OBC — 5k p (0). 
Well below v = 2/i the figure is a bit noisy because the Fourier 
transform in eq. \MiS\ gets its main contribution from a very 
narrow region in space where the solution of the gap equation 
is numerically challenging. 



VI. SUMMARY 

We examined the momentum distribution and momen- 
tum correlations of molecules formed by a Feshbach reso- 
nance or by photoassociation from a quantum degenerate 
atomic gas. Our study elucidated the effect of the atomic 
trapping potential as well as the strength of the atom- 
molecule coupling relative to the characteristic energies 
of the atoms on the molecular momentum distribution. 

Molecules produced from an atomic BEC show a rather 
narrow momentum distribution that is comparable to the 
zero-point momentum width of the atomic BEC from 
which they are formed. In the case of a narrow reso- 
nance, energy conservation limits the molecules to only 
a narrow energy range. The molecule production is a 
collective effect with contributions from all atom pairs 
adding up constructively, as indicated by the quadratic 
scaling of the number of molecules with the number of 
atoms. Each mode of the resulting molecular field is to 
a very good approximation coherent (up to terms of or- 
der 0(1/N)). The effects of noise, both due to finite 
temperatures and to vacuum fluctuations, are of relative 
order 0(1/N). They slightly increase the g^ and cause 
the molecular field in each momentum state to be only 
partially coherent. 

In contrast, the momentum distribution of molecules 
formed from a normal Fermi gas is much broader with 
a typical width given by the Fermi momentum of the 
initial atomic cloud. The molecule production is not col- 
lective as the number of molecules only scales like the 
number of atoms rather than the square. In this case, 
the second-order correlations of the molecules exhibit 
super-Poissonian fluctuations, and it was argued that the 
molecules are well characterized by a thermal field. 
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The case where molecules are produced from paired 
atoms in a BCS-like state shares many properties with 
the BEC case: The molecule formation rate is collective, 
their momentum distribution is very narrow in compar- 
ison to the normal Fermi gas, and the molecular field 
is essentially coherent. The non-collective contribution 
from unpaired atoms has a momentum distribution very 
similar to that of the quasiparticle fluctuations in the 
BEC case. 

In a future publication we will use the Bogoliubov-de- 
Gennes equations to describe the BCS type state which 
is again probed with a molecular field. Thus we will be 
able to go beyond the local density approximation and 
we can study the BEC-BCS cross-over regime. Also, this 
will allow us to study the validity of the local density 
approximation more carefully. 

This work is supported in part by the US Office of 
Naval Research, by the National Science Foundation, by 
the US Army Research Office, and by the National Aero- 
nautics and Space Administration. 

APPENDIX A: CALCULATION OF n(p) FOR 
g > E kin FOR A BEC 

In this appendix we show how the decomposition of the 
atomic field operator into condensed and non-condensed 
fraction together with the local density approximation 
for the non-condensed part can be used to calculate the 
momentum distribution of the molecules for the case of 
an initial BEC of atoms. More details on the properties 
of fluctuations at non-zero temperatures and the use of 
the local density approximation for their description can 



be found in references [U IH, Hi| . 

The starting point of our calculation is eq. © for the 
momentum distribution of the molecules. It reduces the 
problem to evaluating an integral over expectation values 
in the undisturbed atomic field at t = 0. The relevant 
expectation values are most transparently calculated by 
partitioning the quantization volume V into smaller vol- 
umes Vj that are larger than the coherence length £. Par- 
titioning the atomic field operator accordingly, we get 

fa) = J2 ^i( x )Xo(x)c + J2 ^ioc(x), (Al) 

3 3 

where we have introduced functions Aj(x) that are 
one in volume Vj and zero otherwise, and (5^(i) — 
Aj(x)S4>(x). The partitioning of the field operator is il- 
lustrated in fig. [S] We go over to the Fourier transform 

c p = / ^ e -*^(x)=*o(p)2+x; vlN^ ( A2 ) 

where Sc^ (p) is the annihilation operator for a fluctua- 
tion of momentum p in volume Vj. In the local density 
approximation, fluctuations in different cells Vj are un- 
corrected and we have 

(4f £cg> ee 6^,^(5$ Si®). (A3) 

Inserting in eq. © , keeping only terms of first order in 
the fluctuations and making use of relation l|A3|l we find 



n BE c,b(p, t) = {gtfN{N ~l)Y, xS(p/2 - ki)xS(p/2 + k 1 )x (p/2 - k 2 ) Xo (p/2 + k 2 ) 

ki,k 2 

+ %t) 2 A^ ^|xo( P - k)| 2 (c^). (A4) 

I 



The coherent term is readily brought to the form given 
in eq. (JZJ). In the incoherent term we notice that |xo(p)| 2 
is a much narrower function than (cj^c^), the former 
having a typical width of ~ 1/Rtf while the latter has a 
typical width of ~ l/£. Hence, to a good approximation, 
(^k ) can b e treated as a constant for the momentum 



range for which |xo(p)| 2 is nonzero and we obtain 
4( 5 t) 2 iV^||xo(p-k)| 2 <c^c«) 

kj 

«4( ff t) 2 7V^^(c^c«)^|xo(p-k)| 2 

3 k 

= ^) 2 ^7W), (A5) 

3 

where we have made use of the normalization of Xo( x ) 
in the last step. Since £ <C Rtf the volumes Vj can 
be made small and the sum can be approximated by an 
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FIG. 8: Illustration of the partitioning of quantization vol- 
ume and atomic field operator. Also indicated is the function 



integral, leading to eq. JTJl. 



APPENDIX B: CALCULATION OF n(p) FOR 
g < E km FOR A BEC 

The calculation goes along similar lines as in the broad 
resonance case but the evaluation of the integrals is more 
complicated. Repeating the calculation of the broad res- 
onance case that lead to eq. l|A4|l with a p now replaced 
according to eq. (|18fl we find, again to first order in the 
fluctuations 



™bec,„(p, t) = N(N - 1)1 X] A ( E p - e P/2+k - e P /2-k, <)xo(p/2 + k)xo(p/2 - k) 



Vj 



AN J2y \ a ^ e p - £ p/2+k - c P /2-k,*)| Ixo(p - k) 



k,3 



|2 <4 J)t 4 J) > 



^coh 

(p,t) + 

^incoh (p,t). (Bl) 



i 

Let us first consider the coherent part n co h(p, t). Going making the substitution u> — k 2 /M and introducing polar 
over from the summation to an integral in the usual way, coordinates we find 

I 



^coh 



dd sin d 



(P'*) = / <b>A(y-w,t) / 

7r "' z Jo Jo 

x V^Xo (^p 2 /4 + Mu+\p\VMv cos -d^ \o (^Jp 2 /^ + Mlj - \ P \VmZ) cosd^j (B2) 



In the limit t — ► oo, A(u — to) becomes |43 |. 



lim A(f — u,t) = g [ n8{v — u>) — iP 



(B3) 



where, as usual, P means that the integral has to be taken 
in the sense of the Cauchy-principal value. The real part 



J 



r 



can be evaluated by making use of the (5-function and for 
the imaginary part we have to rely on numerical methods 
to calculate the principal value integral. 



Making similar manipulations of the sums over mo- 
menta for the incoherent part leads to 



n inC oh(p, t) = 4^v-^2- V i J M sin * 

<W^|A(i/-o;,i)| 2 |xo ^ p 2 /4 + Mlu - \p\VM^ costij ('(c^c^L-). (B4) 
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Using the delta function to perform the integral in the limit as t goes to inifinity 

we arrive at eq. (|20|l . 
lim |A(z/ - w)| 2 = nq 2 t5(v-u). (B5) 

t — >oo 



S. Inouye et al., Nature (London) 392, 151 (1998). 
E. Timmermans, P. Tommasini, M. Hussein, and A. Ker- [25 
man, Phys. Rep. 315, 199 (1999). 

P. O. Fedichev, Yu Kagan, G. V. Shlyapnikov, and J. T. 
M. Walraven, Phys. Rev. Lett. 77, 2913 (1996). [26 
M. Theis, G. Thalhammer, K. Winkler, M. Hellwig, 

G. Ruff, R. Grimm, and J. H. Denschlag, Phys. Rev. 
Lett. 93, 123001 (2004). [27 
R. Wynar et al., Science 287, 1016 (2000). 
E. A. Donley, N. R. Claussen, S. T. Thompson, and C. E. 
Wieman, Nature (London) 417, 529 (2002). [28 
S. Diirr, T. Volz, A. Marte, and G. Rempe, Phys. Rev. 

Lett. 92, 020406 (2004). [29 
M. Greiner, C. A. Regal, and D. S. Jin, Nature (London) 
426, 537 (2003). [30 
M. W. Zwierlein et al., Phys. Rev. Lett. 91, 250401 [31 
(2003). [32 
S. Jochim et al., Science 301, 2101 (2003). 

M. Greiner, O. Mandel, T. Esslinger, T. W. Hansen, and [33 
I. Bloch, Nature 415, 39 (2002). 
M. Greiner, O. Mandel, T. W. Hansen, and I. Bloch, [34 
Nature 419, 51 (2002). [35 

D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and 

P. Zoller, Phys. Rev. Lett. 81, 3108 (1998). [36 

H. T. C. Stoof (International school of physics 'Enrico 
Fermi', Varenna, 1998). [37 

E. Timmermans, Physica Scripta T110, 302 (2004). 
E. Timmermans, K. Furuya, P. W. Milonni, and A. K. [38 
Kerman, Phys. Lett. A 285, 288 (2001). 

C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. [39 
92, 040403 (2004). 
M. Bartenstein et al., Phys. Rev. Lett. 92, 120401 (2004). [40 
M. Zwierlein et al., Phys. Rev. Lett. 92, 120403 (2004). 
H. T. C. Stoof, M. Houbiers, C. A. Sackett, and R. G. [41 
Hulet, Phys. Rev. Lett. 76, 10 (1996). [42 
G. Bruun, Y. Castin, R. Dum, and K. Burnett, Eur. 
Phys. J. D 7, 433 (1999). [43] 
C. Chin, M. Bartenstein, A. Altmeyer, S. Riedl, 
S. Jochim, J. Hecker Denschlag, and R. Grimm, Science 
305, 1128 (2004). 

J. Kinast, S. L. Hemmer, M. E. Gehm, A. Turlapov, and 
J. E. Thomas, Phys. Rev. Lett. 92, 150402 (2004). 
M. Bartenstein, A. Altmeyer, S. Riedl, S. Jochim, 
C. Chin, J. H. Denschlag, and R. Grimm, Phys. Rev. 



Lett. 92, 203201 (2004). 

E. A. Burt, R. W. Christ, C. J. Myatt, M. J. Holland, 
E. A. Cornell, and C. E. Wieman, Phys. Rev. Lett. 79, 
337 (1997). 

D. Hellweg, L. Cacciapuoti, M. Kottke, T. Schulte, 
K. Sengstock, W. Ertmer, and J. J. Arlt, Phys. Rev. 
Lett. 91, 010406 (2003). 

L. Cacciapuoti, D. Hellweg, M. Kottke, T. Schulte, 
W. Ertmer, J. J. Arlt, K. Sengstock, and L. Santos, Phys. 
Rev. A 68, 053612 (2003). 

E. Altman, E. Demler, and M. D. Lukin, Phys. Rev. A 
70, 013603 (2004). 

R. Bach and K. Rzazewski, Phys. Rev. Lett. 92, 200401 
(2004). 

R. J. Glauber, Phys. Rev. 130, 2529 (1963). 

R. J. Glauber, Phys. Rev. 131, 2766 (1963). 

M. Naraschewski and R. J. Glauber, Phys. Rev. A 59, 

4595 (1999). 

K. E. Cahill and R. J. Glauber, Phys. Rev. A 59, 1538 
(1999). 

D. Meiser and P. Meystre (2004), cond-mat/0410349. 
M. Holland, S. J. J. M. F. Kokkelmans, M. L. Chiofalo, 
and R. Walser, Phys. Rev. Lett. 87, 120406 (2001). 
M. L. Chiofalo, S. J. J. M. F. Kokkelmans, J. N. Milstein, 
and M. J. Holland, Phys. Rev. Lett. 88, 090402 (2002). 
D. A. W. Hutchinson, E. Zaremba, and A. Griffin, Phys. 
Rev. Lett. 78, 1842 (1997). 

J. Reidl, A. Csordas, R. Graham, and P. Szepfalusy, 
Phys. Rev. A 59, 3816 (1999). 

T. Bergeman, D. L. Feder, N. L. Balazs, and B. I. Schnei- 
der, Phys. Rev. A 61, 063605 (2000). 

D. A. Butts and D. S. Rokhsar, Phys. Rev. A 55, 4346 
(1997). 

M. Houbiers et al., Phys. Rev. A 56, 4864 (1997). 

E. Merzbacher, Quantum Mechanics (John Wiley & 
Sons, New York, 1998), 3rd ed. 

We denote the function in Eq. (I19H by A because it has 
properties similar to the usual (5-function. Confusion with 
the local gap parameter introduced below, which we also 
denote by A, cannot arise because the first always has 
energies as its argument while the latter has positions as 
its argument. 



